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SYSTEM AND CIRCUITRY FOR MEASURING ECHOES ON 

SUBSCRIBER LOOPS 

Related Applications 

The present application claims the benefit of U.S. Provisional Application 

5 Number 60/157,094 filed on September 30, 1999, entitled "Method and Circuitry for 
Measuring Weak Echoes on Subscriber Loops". The present application is related to 
US Application No. 09/587,459 filed on June 5, 2000 and entitled "Single Ended 
Method and System for Determining Subscriber Loop Make-Up" and which is 
assigned to the assignee of the present invention. 

io Field Of The Invention 

This invention is related to determining the make-up of copper loops in the 
public switched telephone network and specifically to determining the make-up of 
these loops by processing echo signals generated at discontinuities on the loop. 

Background 

15 The mainstay of the telephone company local network is the local subscriber 

loop. The great majority of residential customers, and many business customers, are 
served by metallic twisted pair cables connected from a local switch in the central 
office to the subscriber's telephones. When customers request service, request a 
change in service, or drop service, these facilities must be appropriately connected or 

20 arranged in the field, referred to as the "outside plant," and telephone companies 
have specially trained craft dedicated full time to this task. Obviously a company 
needs to have an understanding of its subscriber loops including where they are 
connected and the location of the flexibility points such as junction boxes, etc. These 
records historically were kept on paper, called "plats," and more recently are 

25 manually entered into a computer database. However, even when entered into a 
database there are still problems associated with keeping the records accurate and 
up-to-date. 
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Having accurate records of the loop plant is critically important to many 
aspects of a telephone company's business. In addition to the need for accurate 
records to provide traditional voice services, there will be a need for even more 
accurate and detailed records in order to deploy a whole new class of "xDSL" based 
5 services, including those based on integrated services digital network (ISDN), high- 
rate digital subscriber line (HDSL), asymmetrical digital subscriber lines (ADSL) and 
very high rate digital subscriber lines (VDSL) technology. These technologies are 
engineered to operate over a class of subscriber loops, such as nonloaded loops (18 
kft) or Carrier Serving Area (CSA) loops (9 to 12 kft). In fact, the need to be able to 

10 "qualify" a loop for provision of one of these technologies is becoming critical, as 
these technologies emerge and deployment begins. The ability to easily and 
accurately qualify loops will allow telephone companies to offer a whole range of new 
services; problems and high expenses associated with qualifying loops can 
potentially inhibit deployment and/or lower or forego associated new revenues. 

15 Unscreened multipair cables in the existing subscriber loop network constitute the 
main access connection of telephone users to the telephone network. Recently, the 
demand for new services such as data, image and video has increased 
tremendously, and telephone companies plan to deliver broadband ISDN services via 
fiber optic local loops. However, the deployment of fiber optic cables in the access 

20 plant will require at least twenty years, so that, in the meantime, it is extremely 
important to fully exploit the existing copper cable plant. 

Although there are many different digital subscriber line services, for example, 
ISDN basic access, HDSL, ADSL, VDSL, and Synchronous DSL (SDSL), these 
services are not always available to every customer since copper lines present more 

25 problems than expected. In fact, the cable length and the presence of load coils and 
bridged taps may affect the performance of DSL services. Unfortunately, loop 
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records are unreliable and often don't match the actual loop configuration, so existing 
databases cannot be fully exploited. 

Loop prequalification is an important issue not only because it can help an 
economic deployment of DSL services, but also because it can help telephone 
5 companies in updating and correcting their loop-plant records. From this point of 
view, the feasibility of accurate loop make-up identification has much higher 
economic value than simple DSL qualification. 

One way to obtain accurate loop records is to manually examine the existing 
records and update them if they are missing or inaccurate. This technique is 

10 expensive and time consuming. Furthermore, new technologies such as xDSL 
require additional information that was previously not kept for voice services, so there 
is the potential that new information needs to be added to existing loop records. Test 
set manufacturers offer measurement devices that can greatly facilitate this process, 
but typically they require a remote craft dispatch. 

15 Another way to obtain accurate loop records is by performing a loop 

prequalification test. There are essentially two ways of carrying out a loop 
prequalification test: double ended or single-ended measurements. Double-ended 
measurements allow us to easily estimate the impulse response of a loop by using 
properly designed probing sequences. Double-ended testing, however, requires 

20 equipment at both ends of the loop. Specifically, in addition to equipment at the 
Central Office (CO) or near end of the loop, double ended testing involves either the 
presence of a test device at the far end of the loop (Smart Jack or MTU), or 
dispatching a technician to the subscriber's location (SL) to install a modem that 
communicates with the reference modem in the CO. An exemplary double ended 

25 system and method that extrapolates voice band information to determine DSL 
serv i ce capability for a subscriber loop is described in Lechleider, et. al., US Patent 
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6,091,713, entitled "Method and System for Estimating the Ability of a Subscriber 
Loop to Support Broadband Services" (which is assigned to the assignee of the 
present invention). In addition, craft persons may activate softward located at remote 
sites. 

5 In contrast, single ended tests are less expensive and less time consuming 

than double-ended tests. Single-ended testing requires test-equipment only at one 
location, e.g., the CO. In fact, there is no need to dispatch a technician and the CO 
can perform all the tests in a batch mode, exploiting the metallic access with full- 
splitting capability on the customer's line. An example of such a single ended test 

10 system is the "MLT" (Mechanized Loop Testing) product that is included as part of 
the widely deployed automated loop testing system originally developed by the Bell 
System. The MLT system utilizes a metallic test bus and full-splitting metallic access 
relays on line card electronics. By this means, a given subscriber loop can be taken 
out of service and routed, metallically, to a centralized test head, where single-ended 

15 measurements can be made on the customer's loop. The test head runs through a 
battery of tests aimed at maintaining and diagnosing the customer's narrowband 
(4kHz) voice service, e.g., looking for valid termination signatures via application of 
DC and AC voltages. This system is highly mechanized, highly efficient, and almost 
universally deployed. In addition, the MLT system is linked to a Line Monitoring 

20 Operating System (LMOS) thereby providing a means to access and update loop 
records which are useful in responding to customer service requests or complaints. 
However, because this system exclusively focuses on narrowband voice services, 
the system misses important loop make-up features that will be deleterious to 
supporting broadband services via DSL technologies. 

25 Another well known single-ended measurement technique relies on the 

observation of echoes that are produced by medium discontinuities to fully 
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characterize the link. Specifically, these single ended measurements typically rely on 
time domain reflectometry (TDR). TDR measurements are analogous to radar 
measurements in terms of the physical principles at work. TDR test systems transmit 
pulses of energy down the metallic cable being investigated and once these pulses 
5 encounter a discontinuity on the cable a portion of the transmitted energy is reflected 
or echoed back to a receiver on the test system. The elapsed time of arrival of the 
echo pulse determines its location, while the shape and polarity of the echo pulse(s) 
provides a signature identifying the type of discontinuity that caused the reflection or 
echo. Basically, if the reflecting discontinuity causes an increase in impedance, the 

10 echo pulse's polarity is positive; if the reflecting discontinuity causes a decrease in 
impedance, the echo pulse's polarity is negative. A bridged tap, for example, 
produces a negative echo at the location of the tap and a positive echo at the end of 
the bridged tap. Accordingly, a trained craftsperson is able to determine the type of 
fault based on the shape, polarity, sequence of pulses. 

15 Nevertheless, TDR methods (or, in general, single ended measurements that 

rely on echo pulse signatures) are inaccurate and provide ambiguous results that 
even the most skilled craftsperson cannot interpret. Because the arrival of the 
echoes is dependent on the location of the discontinuities (or faults) one echo can be 
masked by another echo if the echoes overlap. In addition, prior art TDR methods do 

20 not take into account, more specifically, are unable to separate, the effects of 
spurious echoes, i.e., echoes generated by a portion of the probing pulse that is 
reflected from a discontinuity, from real echoes, i.e., echoes generated from the 
probing pulse being reflected by a discontinuity. Although spurious echoes will be 
more attenuated than real echoes, they are added to the real echoes causing the real 

25 echo signals between to be distorted. Accordingly, spurious echoes enhance the 
ambiguity inherent in TDR measurements because the shape of the echo is used to 
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interpret the type of fault that caused the echo. In other words, a craft person 
interpreting a TDR measurement analyzes a distorted trace that does not distinguish 
spurious echo distortion. More importantly, the effects of spurious echoes on the 
pulse shape cannot be interpreted via human visual inspection. Further, the effective 
5 range of today's commercial TDRs is quite limited since it is impossible to see 
echoes from discontinuities located more than a few kilofeet away from the detection 
point. Finally, we are unaware of any commercial TDR having the capability to detect 
gauge changes. 

In Galli, et. al. US Application No. 09/587,459 (hereinafter Galli) a method and 
10 system for determining loop makeup based on the echo signatures caused by 
discontinuities as a pulse traverses a loop is disclosed; note that Galli is also a co- 
inventor on the present invention. Although the Galli method is able to determine 
loop makeup more accurately and overcomes the prior art problems highlighted 
above, the method does have some shortcomings. First, the method works well only 
15 where the loop is less than approximately 8,000 feet (8 kft). Once the length of the 
loop increased significantly beyond 9 kft the method is not able to identify loop 
discontinuities with the same success because of noise enhancement due to the use 
of the reciprocal of the insertion loss. Second, Galli's method may not achieve 
unambiguous loop make-up identification if the topology of the loop under test does 
20 not belong to the set of "well-behaved" loops, i.e., loops that follow the 
recommended design rules. 

Of utility then would be a method and system that overcomes the 
shortcomings of the prior art, generally, and more specifically the shortcomings of 
commercial TDRs and Galli. 

25 Summary 

In accordance with the present invention the limitations and shortcomings of 
the prior art are overcome by enhancing the range and resolution of commercial 
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TDRs and by allowing accurate unambiguous determination of a subscriber loop 
make-up. 

In particular, in an aspect of the present invention is a method and circuitry for 
enhancing the performance of commercial TDRs. In more detail, the inductive effect 
5 of a subscriber loop is taken into account when processing echo signals that are the 
result of probing the loop with pulses. In accordance with this aspect of our invention 
we remove a slowly decaying signal caused by the inductive effect of the loop from 
the echo signals. Accordingly, the echo signals are no longer masked by the slowly 
decaying signal thereby increasing the accuracy and range of a measurement 

10 system built in accordance with the present invention. This particular aspect of the 
present invention will be useful in TDR measurements by increasing their range and 
by allowing gauge change detection. The effective range of a TDR designed in 
accordance with this aspect of the present invention is dependent only on the energy 
of the probing pulse, i.e., in principle the range is unbounded. In accordance with 

15 this aspect of our invention, differential probing circuitry that improves the 
performance of TDRs is disclosed. The probing circuitry improves the accuracy of a 
TDR by rejecting the deleterious effects of common mode propagation. 

In another aspect of the present invention, we provide a method for 
unambiguous and precise loop make-up identification. In accordance with this 

20 aspect of the present invention the input impedance of the loop as a function of 
frequency is used in the process of identifying discontinuities and other features 
represented by the echo signals. This is accomplished by first calculating the input 
impedance of the loop as a function of frequency. The input impedance of the loop is 
then convolved, in the frequency domain, with the Fourier transform of the probing 

25 signal. Finally, a simulated waveform of the discontinuity in the time domain is 
obtained by inverse Fourier transforming the result of the convolution. This simulated 
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waveform is then compared to the actual echo signal caused by the discontinuity. If 
the comparison yields an acceptable match, e.g., within a predetermined error 
margin, then the discontinuity is identified and the signal corresponding to that 
discontinuity is removed by subtracting the simulated waveform from the acquired 
5 data. This is done for each discontinuity encountered until the last discontinuity is 
identified. In accordance with this aspect of the invention all the shortcomings of 
Galli are overcome. 

By employing all the foregoing aspects of the present invention, it is possible 
to completely determine the make-up of a loop of any length and any topology. 

10 Brief Description of the Drawings 

FIG. 1 depicts a subscriber loop as a two port network; 

FIG. 2 is the equivalent circuit of FIG. 1 with the loop portion replaced by its 
Thevenin equivalent circuit; 

FIG. 3 depicts the actual voltage appearing across an unterminated, 8 kft 
15 long, 26 AWG loop when a square wave is applied to the loop; 

FIG. 4 depicts the signal of FIG. 3 with particular emphasis on the slowly 
decaying portion of the signal; 

FIG. 5 depicts the echo generated by an unterminated, 8 kft long, 26 AWG 
loop without taken into account a loop's inductive behavior; 
20 FIG. 6 depicts the voltage across an unterminated, 15 Kft long, 26 AWG cable 

when the injected signal is a square wave of 5 microseconds of width and of 1 Volt 
(over 1 00 ohms) of amplitude; 

FIG. 7 depicts echo generated by an unterminated, 15 Kft long, 26 AWG 
cable without taking into account the inductive behavior of the loop; 
25 FIG. 8 depicts the responses of a 8 Kft, 26 AWG unterminated loop for 

different probing signals: square pulse (solid) and half-sine (dashed); 



-8- 



WO 01/24482 



PCT/US00/27052 



FIG. 9 depicts the signals of FIG. 8 with particular emphasis on the slowly 
decaying portion of the signals; 

FIG. 10 depicts the equivalent circuit of FIG.1 used for the computation of the 
slowly decaying waveform superimposed to the echoes; 
5 FIG. 1 1 A illustrates an exemplary loop having a gauge change; 

FIG. 1 1 B depicts the echo response of the loop of FIG. 1 1 A without 
compensating for the inductive behavior of the loop; 

FIG. 1 1 C depicts the echo response of the loop of FIG. 1 1 A after 
compensating for the inductive behavior of the loop; 
10 FIG. 12A depicts the slowly decaying waveform caused by different cable 

gauges over a 50 microsecond time scale; 

FIG. 12B depicts the waveform of FIG. 12A over a 10 microsecond time 

scale; 

FIG. 13 depicts the absolute value of the input impedance versus frequency 
15 for an infinitely long AWG 26 cable; 

FIG. 14 depicts the absolute value of the input impedance versus frequency 
for three different unterminated loops; 

FIG 15 depicts the absolute value of the input impedance versus frequency 
for two loops having bridged taps; 
20 FIG. 16 depicts the absolute value of the input impedance versus frequency 

for two loops having gauge changes; 

FIG. 17A is a high level flow chart that depicts the method steps of the 
present invention; 

FIG. 1 7B is a specific embodiment of the flow chart of FIG. 1 7A; 
25 FIG. 18 depicts the cross-section of a multi-pair twisted cable pair; 
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FIG. 1 9 depicts an equivalent circuit that describes the interaction of the 
propagating modes in the cable of FIG. 18; 

FIG. 20 depicts a prior art application of conventional TDR methods to 
measure one pair of a multi-pair cable; 
5 FIG. 21 illustrates a prior art improvement of the circuit of FIG. 20; 

FIG. 22 illustrates a block diagram of a system that performs broadband 
differential time domain measurements in accordance with an embodiment of the 
present invention; 

FIG. 23 illustrates an implementation of the broadband differential time- 
10 domain sampling head in accordance with an embodiment of the present invention; 

FIG. 24A depicts data acquired on a loop using the measurement setup of 
FIG. 22 in accordance with an aspect the present invention; 

FIG. 24B depicts data acquired on the same loop of FIG. 24A using prior art 
circuitry; 

15 FIG. 25 depicts data acquired, drft), on a real loop using a 3 \xs probing pulse; 

FIG. 26 depicts data acquired, d 2 (t), on the same loop using a 5 \is probing 

pulse; 

FIG. 27 depicts the real loop that the data on which FIG. 25 and FIG. 26 was 
acquired; 

20 FIG. 28 depicts a plot of d } (t) and hj 0) (t) (j = 1,...,4)\ 

FIG. 29 topology 7^ identified in accordance with the method steps of an 

aspect of the present invention; 

FIG. 30 is a plot of e } (t) = dj(t) - h (0) (t) in accordance with an aspect of the 

present invention; 

25 FIG. 31 depicts the families of hypothesized topologies at step i = 1; 

FIG. 32 depicts the hypothesized topologies hj 1} (t) (j=l,...,6) at step i=l; 
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FIG. 33 is a plot of e } (t) and fhj (]) (t) - h (0) (t)] (j = 1,...,6); 

FIG. 34 depicts additional hypothesized topology Tj n at step i = 1; 

FIG. 35 is a plot of e x {t) and h/ l \t) (j = Al, A2, A3); 

FIG. 36 depict cross-correlation functions between e (1 \t) and [h (l \t) - h i0 \t)]-, 
5 FIG. 37 depicts hypothesized topology 7^ at then end of step i = 1; 

FIG. 38 depicts the cross-correlation function between e {l \t) and [A (1) (0 - 

FIG. 39 depicts identified topology 7^ at then end of step * = 1; 
FIG. 40 is a plot of d x (t) and h (1 \t); 
10 FIG. 41 is a plot of e {2) (t) at the beginning of step i = 2; 

FIG. 42 illustrates possible topologies at the beginning of step i = 2; 
FIG. 43 is a plot of e i2 \t) and [h/ 2 \t) - h°\t)] (j = 1, 2); 

FIG. 44 depict cross-correlation functions between e a \t) and [h (2 \t) - h (1 \t)] )] 
using data slice J 7 (0; 

15 FIG. 45 depict cross-correlation functions between e (2 \t) and [/* (2) (0 - h (l \t)] 

using new data slice d 2 (t)\ 

FIG. 46 depicts a chosen topology at the end of step i = 2; and 
FIG. 47 is a plot of d 2 (t), h {2 \t\ and their difference e°\t). 

Detailed Description 

20 The detailed description given below is divided into three sections. In 

particular, section 1 discusses a model for weak echoes, the importance of the 
information contained in the input impedance of a loop, and the method for 
unambiguous loop make-up identification. Section 2 describes a novel differential 
probing apparatus or circuitry that effectively rejects the deleterious effects of 

25 common mode propagation thereby enabling cleaner measurements of echoes. 
Section 3 goes through the process, step-by-step, of loop identification in 
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accordance with the present invention using the method steps of FIG. 17 and the 
measurement setup of FIG. 22. 

1. Method for Unambiguous Loop Make-UP 

a. Weak Echoes: The Underlying Model 
5 In arriving at a method for determining a subscriber loop makeup it is 

necessary to start with a model. The model ultimately determines the success or 

reliability of the method in that the model is an integral part of the method. In the 

present invention (as in the invention of Galli, et. al., US Patent Application No. 

09/587,459 (hereinafter Galli)) the underlying model is critical to the detection of 

10 discontinuities/faults because the model is used in simulating expected responses 
that are compared to actual responses to identify the fault and ultimately determine 
loop makeup. As such, inaccuracies in the underlying model may be amplified as the 
process or method progresses. 

Accordingly, when measuring weak echoes such as those generated by long 

15 loops, the model given in Galli (who is a co-inventor in the present invention) falls 
short. In fact, the previously developed Galli model does not take into account the 
inductive behavior of the loop nor its behavior as a distributed circuit. If an echo is 
very strong, i.e., high amplitude, this behavior may be neglected whereas, in the case 
of weak echoes, neglecting the loops inductive behavior turns out to be too harsh an 

20 approximation leading to incorrect fault detection and limiting the reach of the 
detection system. 

Specifically, lets consider a circuit 100 as shown in FIG. 1, where the loop is 
represented by a two-port network (2PN) 110. If we apply the Thevenin theorem to 
the circuit 100 in FIG. 1, block 120 may be represented by a single impedance equal 
25 to the input impedance of the 2PN (i.e., the input impedance of the loop): 

, „ N AZj +B 

z '» (/)= sf7F (1) 
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The input impedance of the loop is a complex function of frequency and 
obviously depends on the termination of the loop, the load impedance Z L . In the case 
of an unterminated loop (i.e., infinite load impedance), the relationship in eqn. (1) 
reduces to the following expression: 

5 Z<r* m) (/) = lim AZl + B = - (2) 

The waveform actually entering the loop, as shown in FIG. 2, is V 1 and V 1 is 
tied to the source waveform V s as follows: 

As an example let us consider the case where V s is an ideal square pulse of 1 
10 jus width and 1 Volt (over 100 ohms) amplitude injected into an unterminated 8 Kft 
long 26 American Wire Gauge (AWG) cable. The behavior of V^ versus time is 
plotted in FIG. 3. From FIG. 3 it can be seen that the actual waveform entering the 
loop is not a square waveform anymore. In particular, the signal does not go to zero 
at f=1 jus but drops down to approximately 100 mV and then decays very slowly 
15 towards zero. This phenomenon, due to the inductive behavior of a loop, is negligible 
if the amplitude of the echoes is large/strong enough, but has to be taken into 
account when dealing with very weak echoes (low amplitude echoes), in particular, 
the slowly decaying signal and arriving echo signals caused by loop discontinuities 
overlap masking the arriving echo signals. This can be seen in FIG. 4 where an 
20 exploded view of the waveform of FIG. 3 has been plotted from f=20 jus to f=200 jus. 
The spike, 401 , present at approximately 25 jus represents the echo coming back 
from the unterminated end of the loop and is superimposed on the slowly decaying 
signal. The importance of not accounting for the inductive behavior of the loop may 
be better appreciated by reference to FIG. 5, where the echo 501 received from the 
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end of the same loop used in FIG. 4 is shown without being superimposed on the 
slowly decaying signal. As FIG. 5 shows the echo 501 received from the 
unterminated end of the 8 Kft long, 26 AWG loop is prominent and of a much higher 
amplitude than in FIG. 4. Accordingly, where the echo response is weak the 
5 inductive behavior of the loop will mask real echo responses leading to inaccuracies 
in loop makeup identification. 

In principle, the echo overlapping the decaying signal will be visible if the 
echo signal is "sufficiently" strong, has a sharp rise and if it is not too broad. For this 
reason it may be very difficult to detect echoes coming back from very long loops 

10 since those echoes are very weak and broad. As an example, let us consider the 
response of a 15 Kft long, 26 AWG loop to a square pulse of 5 jus of width and 1 Volt 
(over 100 ohms) of amplitude. The echo coming back from the end of the loop on the 
basis of the model in Galli is shown in FIG. 6 whereas the response on the basis of 
eqn. (3) is plotted in FIG. 7. As can be seen from FIG. 7, the echo generated at the 

15 unterminated end is not visible. 

It is worth pointing out that the presence of the decaying waveform is not due 
to the particular choice of the probing signal. In fact, even if a probing signal different 
from a square pulse is used, the received echoes would still be superimposed on a 
slowly decaying waveform since this is due to the intrinsic inductive behavior of the 

20 loop. As an example, let us compare the square pulse with a half-sine pulse (see 
FIG. 8). The half-sine pulse is typically used in high resolution TDRs and it is 
commonly claimed that it yields to higher echo resolution than the square pulse. The 
mathematical expressions of these two probing signals are: 
Square pulse => s sp (t) = A, for0<t<T 

(3a) 

Half -sine pulses Shs (t) = <j2Asm(2ttf 0 t) 9 for0<t<T (f 0 =l/2r) 
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We simulated the responses of a 8 Kft long, 26 AWG unterminated loop to the 
two pulses in Eqn. (3a) and as FIG. 9 shows there is no substantial difference for 
either type of pulse of FIG. 8. The half-sine pulse is widely used in today's high 
resolution TDRs but it does not produce better results than a normal square pulse. 

5 However, there are practical advantages to using a half-sine pulse instead of a 
square pulse. In fact, a half-sine pulse has more energy at low frequencies than a 
square pulse and this property has a twofold advantage. First, it may be more useful 
to detect gauge changes since the reflection coefficient of a gauge change is 
characterized by a low-pass behavior. Second, injecting low frequency pulses in a 

10 pair that's being probed would cause less crosstalk in adjacent pairs that, at the time 
of probing, may be supporting DSL services. Another advantage of using a half-sine 
pulse is that it is easier, from an implementation point of view, to generate "cleaner" 
high-amplitude half-sine pulses instead of high-amplitude square pulses. However, 
other than the above mentioned practical advantages there is no conceptual 

15 difference between the echo response to a square pulse or to a half-sine pulse. 

Since the presence of the decaying signal is unavoidable, the only way to 
reduce or compensate for its effects is to include its effect when processing the 
echoes that result from reflection at discontinuities. Fortunately, it is possible to 
analytically compute the expression of the slowly decaying waveform. In order to 

20 compute this expression, we have to consider that the signal injected into the loop 
« seeS " a | Q ad equal to the characteristic impedance of the first section of the loop. 
This may not be evident from the model shown in FIGS. 1 and 2 because that model 
describes the loop as a 2PN and neglects the fact that the loop is actually a 
transmission line and not a simple circuit with discrete lumped components. The 

25 input impedance of the loop gives a global description of the response of the whole 
loop to a probing signal. But, if we want to take into account that the loop is actually a 
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transmission line, we have to take into consideration that it is the first section of the 
loop that influences the characteristics of the probing signal entering the loop. In fact, 
the probing signal injected into the line propagates along the loop as if it were 
travelling along an infinitely long loop. Since the input impedance of an infinitely long 
5 loop is equal to its characteristic impedance, the correct model describing the voltage 
across the pairs is the one in FIG. 10. This model remains valid until the travelling 
wave encounters a discontinuity along the line, e.g., a gauge change, a bridged tap, 
or an unterminated end. In fact, the presence of a discontinuity along the line causes 
an abrupt change in the boundaries conditions of the equation describing the 

10 travelling wave and, therefore, a change of its shape. Moreover, the changes caused 
by discontinuities on the probing signal travelling along the loop will always take 
place later in time and should not influence at all the shape of the probing signal 
before it encounters these discontinuities. On the basis of the previous 
considerations, we can express the slowly decaying waveform in the following form 

15 (see FIG. 10): 

VAf) = — ^— V s (f). (4) 

Eqn. (4) is the exact expression of the slowly decaying waveform until the first 
discontinuity has been encountered and its effect (the echo) has arrived back at the 
beginning of the loop. 

20 As previously mentioned, the voltage V^f) given in eqn. (3) is the waveform 

obtained when we consider the loop as a discrete lumped circuit (and, therefore, all 
the discontinuities along the loop are included) whereas the voltage V 0 {f) in eqn. (4) 
is the waveform obtained when we take into account the actual nature of the loop as 
a transmission line (and, therefore, only the presence of the first loop section). On the 

25 basis of the previous considerations, subtracting V Q (f) from !/,(/) should remove the 
slowing decaying waveform due to the inductive behavior of the loop and allow for 
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easier detection of weak and broad echoes. So, the expression, V(f), for an echo 
signal, taking into account the effects of the slowly decaying signal (i.e., the inductive 
behavior of the loop) should be approximately given by the following expression: 

Z ; „ Z„ 



V(f) = V,(f)-V 0 (f) = 



z in +z s z 0 +z s 



v s (f) = 

(5) 



(Zi^JZs .y s(f) 



(z in +z s )(z 0 + z s ) 

5 The subtraction of waveform eqn. (4) from waveform eqn. (3), resulting in 

eqn. (5), is a very useful technique in that it allows us to detect discontinuities located 
far away from the measurement location. Moreover, this technique is also very useful 
for the detection of gauge changes which are not detectable unless the slowly 
decaying waveform is removed as by eqn. (5). As an example, let us consider the 

10 echo response of the loop in FIG. 1 1 A that contains a gauge change 1 101 . The echo 
response without accounting for the decaying signal (using eqn. (3)) is shown in FIG. 
11 B. Even though the discontinuity is very near (just 1000 feet away from the 
measurement equipment), only the echo 1110 given by the unterminated end of the 
loop is clearly visible. However, using eqn. (5) we obtain the plot in FIG. 11C which 

15 clearly shows the presence of another echo 1113 that is due to the gauge change. 

In principle, the computation of the slowly decaying waveform in eqn. (4) 
requires the knowledge of the type of gauge of the first loop section. The gauge of 
the first loop section is constituted by the cable coming out of a Central Office (CO) 
and it is reasonable to assume that this information is known. However, if this 

20 information were not available a priori it is still possible to identify the kind of gauge of 
the first loop section by trial and comparison. FIG. 12 shows the behavior of eqn. (4) 
for different kinds of gauges and for a 1 Volt and 100 ns square pulse probing signal. 
As it can be seen from FIG. 12, different kinds of gauges produce different slowly 
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decaying waveforms. This suggests that it is convenient to probe the loop with very 
short pulses and, then, compare the observed signal with the signals given by all the 
possible gauges. In fact, looking at FIG. 12A we can see that the waveforms tend to 
be very similar after the first 20 microseconds, whereas they are quite different in the 
5 first 15 microseconds (see FIG. 12B). This property can be exploited by computing 
the difference between the observed slowly decaying waveform and the simulated 
ones and, then, determining the kind of gauge that yields to the smallest difference 
signal. 

Those or ordinary skill in the art will note that the waveform or signal 

10 corresponding to V 0 {f) in eqn. (4) can be obtained by probing a long loop, e.g., 20 kft 
- 30 kft, of characteristic Z 0 . As a practical matter, the long loop can be probed using 
a conventional or prior art TDR; that is, V 0 (f) in eqn. (4) can be obtained empirically. 
The result of the probing can be stored in memory and subtracted from the waveform 
corresponding to the loop under the test, i.e., the acquired data, which would be 

15 obtained by probing the loop under test. Subtracting this empirically waveform from 
acquired data would be equivalent to subtracting a simulated waveform from the 
acquired data as indicated in eqn. 5. 

It is worth pointing out that this technique has never been used in commercial 
TDRs. Vendors claim that today's TDRs have ranges up to 50,000 feet and a 

20 resolution under one foot. However, the reality is that neither of these claims are 
achievable without first removing the effect of the slowly decaying signal on the 
measurements. From this point of view, the previously described technique could be 
extremely useful to improve the performance of a TDR. in particular, removing the 
effect of the slowly decaying signal and using the circuitry described in section 2 

25 yields superior performance over prior art TDRs. 
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b. The Information Contained in the Input Impedance of a Loop 

Removing the slowly decaying signal, in effect compensating for the inductive 
5 behavior of the loop, as described above in section 1a allows us to obtain much 
cleaner responses than those available using commercial TDRs. In particular, this 
technique has allowed us to detect gauge changes and discontinuities located far 
away from the detector. Another important result of section 1a is that the input 
impedance of a loop can be viewed as the frequency domain description of the echo 

10 observation process. In this section we will describe how the model of section 1a, 
i.e., the frequency domain representation of the loop input impedance that includes 
the inductive behavior of the loop, can be used to identify the makeup of a loop and 
detect discontinuities located far away, including gauge changes. 

Since a frequency-domain reflection measurement is the composite of all the 

15 signals reflected by the discontinuities over the measured frequency range, the input 
impedance of a loop contains information on both the real and the spurious echoes 
that would be generated by probing the medium with a signal. In fact, in our work we 
have found that different loops or different loop discontinuities provide recognizable 
signatures in the loop input impedance. The behavior of the absolute value of the 

20 input impedance of a few exemplary loops is plotted versus frequency in FIGS. 13 - 
16. In particular, the input impedance of an infinitely long 26 AWG cable is shown in 
FIG. 13. As FIG. 13 shows, as the frequency increases, \Z in \ decreases monotonically 
to the asymptotical value of 100 ohms. 

On the other hand, if a discontinuity is present along the line, the behavior of 

25 \Z in \ no longer monotonically decreases and, moreover, different discontinuities yield 
to different behaviors of the input impedance. Specifically, let us consider the case of 
an unterminated loop as is shown in FIG. 14. In this case, \Z in \ initially starts to 
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decrease monotonically and, then, exhibits a damped oscillatory behavior starting at 
a few tenths of khz. The amplitude and the frequency of these oscillations are tied to 
the length of the loop: for short loops the amplitude is higher and the frequency is 
lower and viceversa for longer loops. In FIG. 14 we show the behavior of \Z in \ versus 
5 frequency for three type of unterminated loops: a 1 kft long AWG 24 unterminated 
loop 1401, a 10,5 kft long AWG 24 unterminated loop 1411, and a 15 kft long AWG 
26 unterminated loop 1421. 

Turning now to FIG. 15, there is shown \Z in \ for a loop having a bridged tap. 
As FIG. 15 shows, \Z in \ initially starts to decrease monotonically and, after a few 

10 tenths of khz, the line exhibits a recurrent pattern: a small peak followed by a bigger 
peak, both with damping amplitude. However, since the frequency and the amplitude 
0 f these peaks depend on the lengths of the loop sections constituting the bridged 
tap, it is very difficult to establish a general pattern for all the possible configurations 
of a bridged tap. Specifically, in FIG. 15 we show the behavior of \Z in \ versus 

15 frequency for two cases: an AWG 24 loop having a 980 ft long bridged tap located 
980 ft on the loop 1501 and an AWG 26 loop having a 1500 ft long bridged tap 
located 6000 ft on the loop 1511. 

Finally, the case of a gauge change is shown in FIG. 16. As shown in FIG. 
16, a loop with a gauge change exhibits a ripple but this ripple is very small and 

20 difficult to detect. This is not a surprise because, among all the kinds of 
discontinuities we considered, the gauge change is indeed the hardest discontinuity 
to detect. 

Based on the foregoing, we note the following features of different loops or 
loop discontinuities. A common property to all the discontinuities is that, when a 
25 discontinuity is near, the oscillatory behavior is much more pronounced than for the 
case of a discontinuity located farther away. This suggests that the "onion peeling" 
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or "de-embedding" method of Galli may also be used to detect/identify discontinuities 
by sequentially removing the echoes based on their position on the loop, starting with 
the discontinuity located closest to the measurement equipment and ending with the 
farthest discontinuity. More importantly, the loop impedance frequency responses of 
5 the exemplary loops considered above indicate that the loop can be considered a 
deterministic or known channel in the frequency domain. On the other hand, the 
behavior of the input impedance of a loop where discontinuities are located far away 
resembles the behavior of an infinitely long loop. This property suggests a dual 
interpretation in the frequency domain of the difficulty of detecting far discontinuities. 

10 In fact, note that in the time domain discontinuities located far from the point of 
detection are very hard to detect because the echoes coming back are very small 
and very broad, whereas, in the frequency domain, discontinuities located far from 
the detection point are difficult to detect because they exhibit almost the same 
behavior of an infinitely long loop, i.e. a loop without discontinuities. 

15 The deterministic nature of the loop had been exploited in the identification 

method proposed by Galli and particularly in the active loop make-up identification 
phase. Galli proposed the use of the reciprocal of the insertion loss of the identified 
echo path in order to limit the broadness of the far echoes. However, as previously 
pointed out that method is very sensitive to noise enhancement and this limits its 

20 application to loops up to 8 - 9 kft long. In fact, in Galli we estimated the noise power 
to be approximately -120 dBm/Hz but the real noise power on a loop (including 
crosstalk and other sources) may be even higher so that the actual range may be 
even lower. In accordance with the present invention, the problem of the noise 
enhancement is avoided because the reciprocal of the insertion loss is not used to 

25 determine the loop makeup. 
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c. Justification for the Proposed Approach 

The waveform observed at the receiver after probing a loop with a signal is 
constituted by an unknown number of echoes, some overlapping some not, some 
5 spurious, some not, that exhibit unknown amplitude, unknown time of arrival and 
unknown shape. This problem is very complicated and has been seldom addressed 
in scientific literature. Due to these analytical difficulties, we took a different approach 
to solving the problem. 

Given the availability of an accurate model of the physical phenomenon, the 
10 most reasonable choice for carrying out the identification process is to apply the 
Maximum Likelihood (ML) principle. The ML method is based on the idea that 
different populations generate different samples and that any given sample is more 
likely to have come from some populations than from others. 

Similarly, our method hypothesizes a set of loop topologies and, on the basis 
15 of the mathematical model, computes the waveforms that should be observed at the 
receiver if the hypothesized topologies were true. The topology corresponding to the 
waveform that best matches the observed signal is chosen as the "most likely" 
topology. An index of the "closeness" between the hypothesized waveform and the 
observed one may be, for example, the Mean Square Error (MSE). In principle, if the 
20 model were exact and no form of noise were present, this procedure is able to 
identify exactly and flawlessly a loop. Obviously, noise is always present and the 
mathematical model, although precise, is not exact. Therefore, in practice, there is 
the possibility of making errors. An assessment of the probability of erroneous 
identification necessarily requires an extensive measurement campaign, especially 
25 field trials were impulse noise, crosstalk and non-ideal situations are often present. 

However, in practice, the ML principle cannot be applied exactly as previously 
stated. In fact, it is very impractical to hypothesize a certain loop topology since the 
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set of all the possible loops is too vast. A loop is constituted of several loop-sections 
made of different gauges and spliced together in a certain number of ways. There are 
only four possible gauge kinds (AWGs 26, 24, 22, 19) and only four main 
discontinuities (gauge change, bridged tap, co-located bridged taps, end of loop). 
5 The number of loop sections constituting a loop is limited to not more than 9 or 10 
sections. However, the location of each discontinuity, i.e., its distance from the CO, is 
a parameter that can assume a set of non-numerable values. This is why it is 
impossible to apply the ML principle previously described to the set of all possible 
topologies. 

10 A way to avoid this problem is to follow a step-by-step approach and apply 

the ML principle at every step. More in detail, discontinuities should be identified one 
at a time starting from the nearest one and ending with the last one. So doing, the ML 
principle will be applied at every step to identify a single discontinuity. In this way, the 
topology of the loop under investigation will be identified one section at a time and 

15 loop sections will be added to the hypothesized topology one at a time as well. 
However, the identification of a discontinuity will completely characterize the topology 
only up to the last identified discontinuity, whereas nothing could be said about 
whatever follows the last identified discontinuity. This implies that the loop section 
following the last identified discontinuity will not be identified in terms of length or 

20 gauge since this could only be accomplished when the discontinuity following that 
section will be observed and identified. This is a big problem since there is no way 
that the next discontinuity can be identified if no hypothesis is made on the loop 
section that precedes that discontinuity. However, a simple way of solving this 
problem is to hypothesize that the loop section that follows the last identified 

25 discontinuity is constituted of an infinitely long loop. This will allow us to identify the 
following discontinuity since the behavior of that loop section will certainly be different 
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from the behavior of a loop section that has no discontinuities at all. Moreover, the 
hypothesis of an infinitely long loop section will also be useful in the elimination of the 
always present slowly decaying signal. 

This step-by-step approach is advantageous for at least two reasons. First, 
5 the nearest discontinuities are the easiest to detect since they generate stronger 
echoes and these echoes are not hidden by other echoes. Secondly, once a 
discontinuity has been identified, its echoes (both real and spurious) can be 
subtracted from the observed waveform so that the two important problems of 
overlapping echoes and of the presence of spurious echoes are simultaneously 
10 solved. 

Interestingly, the ML procedure adopted for our method is similar, in principle, 
to the ML estimation of the sequence of the states of a Markov chain in additive 
noise. As it is well known, a recursive solution to this problem is given by the Viterbi 
algorithm. The Viterbi algorithm computes all the admissible transmitted sequences 

15 and chooses as the transmitted one the one that is "nearest" to the observed signal 
in an Euclidean sense. The Viterbi algorithm is optimal in the sense that the decided 
sequence is the nearest sequence to the received one. From this point of view, it is 
reasonable to say that our method is optimal in the sense that the decided topology 
is the topology that generates the waveform that is nearest to the observed one. It is 

20 worth pointing out that this is not the same thing of stating that the decided topology 
is the nearest topology to the real one. In fact, in order to formally prove this it is 
necessary to introduce a metric in the space of the loop topologies and the definition 
of this metric is not straightforward. 

d. Method Steps for Loop Identification 
25 Turning now to FIG. 17A, there is shown a high level flow chart depicting the 

essential method steps of the present invention. In general, given a set of acquired 
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data from step 1701, our method proceeds as follows we hypothesize a set of loop 
topologies 1702 and, on the basis of the mathematical model, compute the 
waveforms that should be observed at the receiver if the hypothesized topologies 
were true 1703. The topology corresponding to the waveform that best matches the 
5 observed signal is chosen as the "most likely" topology 1 704. 

Turning to FIG. 17B there is shown a specific embodiment of the high level 
flow chart of FIG. 17A. Specifically, at step 1701 the process begins by acquiring 
data. The process for acquiring data essentially comprises repetitively sending 
pulses of varying pulse widths onto the loop and receiving the echo responses 
10 generated by discontinuities on the loop, i.e., the signal of the echo responses 
comprising the acquired data. In our notation the acquired data is represented by 
vector d(f). 

It is necessary to probe the medium with probing signals of different width for 
severa i reasons. The width of the pulse determines our blind zone, i.e., the zone 

15 where discontinuities cannot be detected. In fact, it is impossible to detect echoes as 
long as the probing pulse lasts. This implies that narrow pulses are necessary to 
detect near discontinuities. On the other hand, broad pulses are necessary to detect 
far discontinuities, since there is not enough energy in narrow pulses to reach the 
end of long loops. We have ascertained from our measurements that a square pulse 

20 of different widths is sufficient to identify loops within 18 kft. In our work we have 
found pulses having widths of 500 ns, 1 jus, and 5 jus to be adequate. The 500 ns 
pulse implies that the first discontinuity that can be detected is located at 
approximately 160 ft from the CO. 

Square pulses provided very good results in our measurements. However, 

25 their use may be limited in the field if certain switches are present in the CO. In fact, 
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some switches may exhibit a low-pass filtering behavior so that pulses with narrower 
bandwidth should be used in these situations. 

From the foregoing, during the data acquisition phase 1701 several 
snapshots of the medium will be taken. In general, the snapshots corresponding to 
5 the narrower pulses will be analyzed first and the snapshots corresponding to the 
broader pulses will be analyzed afterwards. 

Once the data is acquired, the next step is data processing 1705. Data 
processing comprises many steps as FIG. 17 shows. We will now describe each of 
those steps in detail. 

10 The first step in processing the data is identifying the gauge of the first loop 

section 1708, i.e., the loop section attached to or immediately following the CO or 
measurement equipment. In accordance with our method identifying the first loop 
section gauge comprises substeps 1709, 1711, and 1712. 

At substep 1709 an infinitely long loop is hypothesized. In so doing, four loop 

15 topologies are hypothesized, one for each gauge kind, i.e., one each for AWG 19, 22, 
24, and 26 cable. In general, any slice can be used as long as the metrics are 
computed on the slowly decaying signal only. The loop impedance response of the 
infinitely long loop is based on the model discussed previously in section 1.a; in fact, 
refer to FIG. 14 for examples of the loop impedance response for infinitely long loops. 

20 At step 1711 the slowly decaying waveform for each of the hypothesized 

infinitely long loops is then calculated, block 1711. The details of calculating the 
slowly decaying waveform of a loop have been fully explored in section 1 .a. 

Once the slowly decaying waveform for each hypothesized infinitely long loop 
is calculated, as in step 1711, the slowly decaying waveform that best matches d(t), 

25 or the data slice, is chosen as the waveform representing the first loop section, step 
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1712. In effect, the gauge of the slowly decaying waveform that best matches the 
data slice is chosen as the gauge of the first loop section. 

With the gauge of the first loop section identified, step 1 708, the method then 
proceeds to step 1713. At 1713 the slowly decaying signal (SDS) or waveform 
5 (SDW) that represents the first loop section is subtracted from d(t). That is, if in 
general h {i) (f) is the simulated waveform corresponding to the identified topology 7^ 
at step /, then e (i) (t) = d(t) - h {l ~ 1) (t), the difference between the acquired or real data 
waveform and the simulated waveform corresponding to the identified topology 7^' 1) . 
At this step in processing, prior to processing the echo corresponding the first 

10 discontinuity, h {l ~ 1) (t) would correspond to the slowly decaying waveform that best 
matches the data slice; specifically, / = 1 , and e (1) (t) = d(t) - hf 0) (t), where h (0) (t) 
represents the chosen SDS at 1712. As such, e {f \f) represents the acquired 
waveform without the slowly decaying waveform. At this stage of the process e {i) (t) 
can be thought as being representative of compensated data, i.e., the data has been 

15 adjusted to compensate for the inductive behavior of the loop. 

With the effects of the decaying waveform, or conversely the effects of the 
inductive behavior of the loop, compensated for, the remaining task essentially 
comprises looking at the echo responses making up the remaining acquired data, 
matching those acquired echo responses to echo responses predicted by the model 

20 of section 1.a, and, once the predicted echo response matches the actual echo 
response within an error margin, identifying the discontinuity that causes the echo. 

Accordingly, processing continues with estimating the time of arrival of the 
first echo, step 1714. The estimation of the time of arrival of an echo is of 
fundamental importance. In fact, since the velocity of propagation of an electrical 

25 signal on a twisted-pair is known (approximately equal to 0.66 c, where c is the 
speed of light in vacuum) the knowledge of the time of arrival of an echo allows us to 
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determine the location of the discontinuity. If r is the time of arrival of the detected 
echo, the location of the corresponding discontinuity may be computed as follows: 

7 1 

I = — V-T , 

2 

where v is the velocity of propagation of the echo and the factor Vi takes into account 

5 the fact that the signal has traveled round-trip from the CO to the discontinuity. In 
more detail, the speed of an electrical signal on a twisted-pair depends on the kind of 
gauge, although these speeds are very close. In general, signals travel faster on 
thicker cable than on thinner ones. So, when estimating the location of a discontinuity 
a precise estimate would involve the calculation of the travelling time of the echo 

10 along each different loop section, i.e., considering different velocities of propagation 
for each loop section. 

As the identification process goes on, it is possible that some small errors in 
the estimation of the location of the first discontinuities add up and cause larger 
errors in the estimation of the location of the following discontinuities. This fact may 

15 not be a serious problem, since errors of a few hundred feet are certainly not critical 
on loops as long as several thousand feet. However, if this error grows too much it 
may cause problems in the choice of the correct topology and jeopardize the correct 
identification of the loop. For this reason, it is important to accurately estimate as 
much as possible the time of arrival. We have found in our work that process is more 

20 accurate if the time of arrival estimation is done at two different moments. First, a 
rough estimate is performed at step 1714 and, afterwards, this estimate is refined at 
step 1740. It is noted that step 1740 is optional because it may be possible to 
accurately estimate the time of arrival at step 1714 using more sophisticated data 
processing techniques. Furthermore, accuracy is a relative concept and we have 

25 found that additional accuracy in the time of arrival did not lead errors in loop make- 
up identity in the vast majority of cases we tested. 

-28- 



WO 01/24482 



PCT/US00/27052 



With the time of arrival estimated we then hypothesize a representative set of 
reasonable topologies that could match the next discontinuity 1715. In our notation 
the set of all hypothesized topologies for the / h echo, or iteration through the method 
steps, is represented by {7} (0 }. There are a finite number of possible topologies that 

5 can be hypothesized, and this number is also small. This suggests that a simple 
exhaustive search through all the possible topologies could be performed without 
requiring a prohibitive computational burden. However, several tens of topologies 
may be possible in some cases such as several co-located bridged taps. In this case, 
it could be useful to organize the topologies in "families" or clusters and chose the 

10 two topologies (sample topologies) that exhibit the most "distant" waveforms. In this 
way only two topologies per cluster would be tested, thus reducing the computational 
burden. Once this preliminary test has been performed, the sample topology that 
best matches the observations will define the family of topologies that most likely 
contains the best topology. At this point we would limit the search for the best 

15 topology within that family of topologies. The set of hypothesized topologies 
depends on the sign of the detected echo and on the previously identified 
discontinuity. 

Obviously, if some a priori knowledge on the statistical distribution of the loop 
sections is known, the search could be performed more efficiently. In fact, the most 
20 recurrent topologies may be hypothesized first so that the time needed for the 
determination of the gauge or for the determination of the kind of discontinuity can be 
reduced. A partial statistical characterization of the loop plant could be obtained by 
analyzing the loop records of the CO under test contained in the LFACS database. 

Also the discontinuities encountered before the discontinuity currently under 
25 identification determine the possible topologies that can be identified. In particular, 
two cases have to be considered: the case when a previous bridged tap has already 
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been thoroughly identified (both in location and length) or if the bridged taps have 
only been located but their length is still unknown. In the following, these two cases 
will be referred to as "closed" or "open" bridged tap cases, respectively. The reason 
why it is important to discriminate between the two cases of open and closed bridged 
5 taps is due to the fact that a bridged always generates a pair of echoes, a negative 
echo followed by a positive echo. However, it is not necessarily true that the positive 
echo immediately follows the negative one. For this reason, once a bridged tap is 
located it has to be taken into account that the next observed echo might be 
generated either by the end of the bridged tap or by some other new discontinuity. 

10 Finally, all the loop sections share the property of being considered 

constituted of infinitely long sections that for now are of unknown length. 

A simulated waveform that represents each hypothesized topology is then 
generated 1717 according to the mathematical model of section 1 .a. The process for 
simulating a waveform for a topology 1717 is given by steps 1718 through 1720. 

15 First, the input impedance of the each hypothesized loop topology is calculated as a 
function of frequency 1718. The calculated input impedance of each hypothesized 
loop topology is then convolved, in the frequency domain, with the Fourier transform 
of the probing signal 1719 as indicated by the expression of Vrff) in eqn. 5. Finally, a 
simulated waveform for each hypothesized loop topology is obtained in the time 

20 domain by inverse Fourier transforming the result of the convolution 1720. As 
previously indicated, by our notation {/)/*(*)} represents the set of all the simulated 
waveforms (/=1, A/ /} ) corresponding to the set of all the possible hypothesized 
topologies {7j (0 } 

We then compare the simulated waveforms of step 1717 to remaining 
25 acquired or observed data signal, e (0 (f)> and choose the simulated waveform that 
best matches e (/) (0, step 1731. There are many ways to determine what constitutes 
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the "best match" or has the smallest error margin. The topology corresponding to the 
waveform that best matches the observed signal is chosen as the "most likely" 
topology. An index of the "closeness" between the hypothesized waveform and the 
observed one may be, for example, the Mean Square Error (MSE). In principle, if the 
5 model were exact and no form of noise were present, this procedure is able to 
identify exactly and flawlessly a loop. Obviously, noise is always present and the 
mathematical model, although precise, is not exact. Therefore, in practice, there is 
the possibility of making errors. In our work we chose the MSE as the metric for 
comparing the real data waveform and the simulated ones. A set of A/° MSEs is 

10 computed between the acquired data waveforms and the simulated data waveforms, 
i.e., e (i) (t)=(d(t) - h (i) (t)), and the ^ differences ({h/ i} (t)} - h (i - l) {t)) between the 
waveforms pertaining to all the discontinuities topologies at step / and the waveform 
corresponding to the discontinuity topology identified at the (/-1)-th step. 

The waveform obtained at step 1731 is then subtracted from the remaining 

15 acquired data 1750, i.e., de-embedding. Thus, in accordance with our we out 
notation the remaining signal would be given by e (i) (f) = d(t) -h {l ~ 1) (f). 

The process then continues by checking for the existence of another echo, 
1760. If the signal level indicates that the existence of another echo the process 
returns to step 1714 otherwise a loop record is created as at step 1764. 

20 Note that, similar to Galli, the basic idea of proceeding step-by-step in the 

identification of the loop discontinuities and of stripping off the effect of near 
discontinuities from the far ones is still valid, but the way these effects are stripped off 
from the observation is different. This "de-embedding" technique is more effective 
than that described in Galli because it would not suffer the problems of noise 

25 enhancement due to the use of the reciprocal of the insertion loss. Moreover, since a 
frequency domain reflection description is the composite of all the signals reflected 
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by the discontinuities over the measured frequency range, the signal that is 
subtracted from the acquired data will also contain all the spurious echoes generated 
up to the last identified discontinuity. This avoids the use of the time domain model 
for spurious echoes where all the spurious echoes have to be generated one by one 

5 and then subtracted from the acquired data order to be removed. 

This method, in contrast to Galli which uses the reciprocal of the insertion 
loss, will not modify the broadness of the received echoes. However, this is not a 
problem because, after having removed all the previous echoes, we are sure that the 
echo pertaining to the next discontinuity will be the first one to appear and will not be 

10 hidden by the previous ones. It is important to point out that, since this method uses 
a frequency domain model, the IFFT of the input impedance will contain both the real 
and the spurious echoes, thus making the de-embedding complete. Finally, in 
accordance with this aspect of our invention and in contrast to Galli we can identify 
loops of any topology. 

15 2, Circuitry for Improving the Detection of Weak Echoes 

The method described in the preceding section 1 improves the reach and 

accuracy of single ended metallic measurements methods by formulating a more 
accurate model of the behavior of the signals that are used to make the 
measurement. Specifically, the method took into account the inductive behavior of 

20 the loop and compensated for its effect on the measurements. In addition, the 
previously described method used the frequency domain representation of the 
discontinuities to overcome noise enhancement in the time domain. During the 
course of our work, we also found that we can further improve performance through 
the use of differential probing. In this section 2 we will now describe our differential 

25 probing circuitry and provide an example of the improvements that may be obtained 
by use of our circuitry when used in conjunction with the method of section 1 . 
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10 



Multi-pair twisted pair cables exhibit more complex behavior than simple 
cables often characterized using time-domain reflectometry. FIG. 18 illustrates the 
cross-section of a multi-pair twisted pair cable 2100 which includes a test pair 
consisting of wires 2102 and 2104, a plurality of inactive pairs 2106 and a ground 
path 2108. It is desirable that the inactive pairs do not affect measurements 
performed on the test pair. Wires 2102 and 2104 have voltages V wl and V w2 relative 
to ground 2108, respectively. The inactive pairs 2106 are assumed to be 
equipotential at voltage V w3 for the purposes of analysis. This system of conductors 

can support three propagating modes (low frequency or Transverse Electric and 
Magnetic (TEM) approximation) which can be described by: 
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The factor 6 describes the shielding produced by the inactive pairs: complete 
15 shielding gives 0-0, while 9 =0.5 for twisted pair cable. It is useful to note that 
B~ l =A T and A~ l =B T . Other test pair combinations can be treated by 
superposition. 

FIG. 19 is an equivalent circuit which describes the interaction of the 
propagating modes in cable 2100 represented by eqns. 6(a) and 6(b). The terminals 
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corresponding to the distal end of the cable are represented by nodes 2202, 2204 
and 2206. V wl and I w] denote the voltage and current on wire 2102 (w1), etc. The 
differential mode I dif represents current confined to wires 2102 (w1) and 2104 (w2) 
and is generally the desired signal. The mode I pr represents current flowing 
5 between the plurality of inactive pairs 2106 (w3) and the test pair of wires 2102 and 
2104. The two modes I dif and I pr are well confined to the cable and consequently 
exhibit low attenuation. The characteristic impedances associated with these modes 
are typically Z^=100 - 900 Q and Z pr =20 - 75 £2, respectively. The net cable 

current / is another common-mode signal associated with both radiated emissions 
10 and noise pick-up, and is highly dependent on cable installation. The characteristic 
impedance of this lossy mode Z cm is variable and not readily characterized. 
Excitation of the inter-pair mode and cable common mode can adversely affect 
measurement of differential mode signals. 

There are six propagating voltages and currents related by Vj + = Z 0 /j + and 

15 V~ =Z 0 1~ , where Z 0 is the diagonal matrix of characteristic impedances. 

Discontinuities in the cable due to imperfections, interconnections or cable 
mismatch can induce undesirable coupling between all three modes. Consider a 
semi-infinite cable which is excited by a purely differential incident traveling-wave 

Il = (l dif ,0,0) r and whose termination is represented by an impedance matrix Z term . 
20 The boundary condition in this case is V cl = Z term I cl . The reflected wave /" may 
contain additional common-mode components I pr and I cm . Expressed in terms of 
the propagating modes using Eq. (6), the current and voltages satisfy I^-I^= BI cl 
and V* + V } ~ = AZ /erffl I cl at the termination. Combining these relations yields: 
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(Z 0 +AZ term A T )I-=(AZ temi A T -Z 0 )ir . (7) 
Standing waves can occur with each of the three fundamental modes and adversely 
affect measurements. Reflections are suppressed when AZ term A T = Z 0 . Using the 

symmetry between A and B noted earlier, this can be expressed as Z (erm = B T Z 0 B, 

5 which is just the impedance looking into an infinitely long cable. 

Another important case is mode coupling in a cable interrupted by co-located 
impedances in series with each wire. This is equivalent to the termination 
Z term =Z se +B T Z 0 B where Z se is a diagonal matrix whose elements are the series 

impedances. Using Eqn. 7, the reflected wave I~ in this case satisfies 

10 (Z 0 + l/2AZ se A T )I" = 1/2A Z se A T 1+ . (8) 

The transmitted wave in the cable beyond the position of Z se is / + = /+ . Eqn. (8) 

has the iterative solution 7" = £ (- 1) (l/2 Z 0 _1 AZ se A T ) 7, + . 

For a single resistance R se in series with wire 1 or wire 2, Eqn. 8 reduces to 



V _ R se 



i+ — —+ — —+ 

2Z dif 8Z pr 8Z cm y 



(9) 



15 which has little dependence on Z cm , and 

i ez 

_£HL = EL , (10) 

■^pr ^cm 

which is independent of R se . 

A cable interrupted by localized shunt conductances between the wires and to 
ground can be described by / cl = / c2 + ^ v a with ^1=^2^ where the subscripts d 
20 and c2 denote the two sections of the divided cable. One can show in this case that 
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a+y 2 BY sh b t z 0 )i: = i+ (id 

and = . A single shunt conductance Y sh from wire 1 to ground produces a 

common-mode current I cm II + dlf ~ Y sh Z dlf /4, which does not benefit from shielding by 

the other pairs. A shunt conductance placed between the conductors does not 
5 degrade CMR. Based on Eqn. 1 1 , the common-mode rejection is: 

hlIL = - Z se , Y *h Z dif Z 3 (12 ) 

i; f "2 Z cm +Z 3 2 z cm +z 3 
FIG. 20 illustrates an application of conventional TDR methods to measure 
one pair of a multi-pair cable. This figure includes loop 2100 and a test system 2301 
consisting of pulse generator 2310 and oscilloscope 2320. The test system has 
10 terminals 2330 and 2332 connected to wires 2102 and 2104 of loop 2100, 
respectively. Pulse generator has single ended output 2312 which is connected 
through an output impedance Z 0 2316 to terminal 2330, and signal return 2314. 

Oscilloscope 2320 has input 2322 connected to terminal 2330, and signal return 
2324. There is a system ground 2302 to which both signal returns and terminal 2332 

15 are connected. System ground 2302 may not coincide with cable ground 2108. With 
this method, an unbalanced signal is applied to wires 2102 and 2104 which excites 
all three cable modes described previously in eqns. 6-12. The oscilloscope also 
responds to signals from all three modes. 

FIG. 21 illustrates an improvement of the prior art that partly alleviates 

20 common mode interference. FIG. 21 includes the circuit of FIG. 20 with the addition 
of transformer 2410 and termination resistor 2420. The primary side of transformer 
2410 is connected between terminals 2330 and 2332. The secondary of transformer 
2410 has balanced outputs 2412 and 2414 and center tap 2416. The balanced 
outputs are connected to loop access terminals 2102 and 2104. Center tap 2416 is 
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connected to system ground 2302 through resistor 2420. Transformer 241 0 converts 
single ended signals appearing at terminal 2330 into a differential signal applied 
between 2102 and 2104. Differential signals are also converted to single ended 
signals for display on the oscilloscope. Resistor 2420 can terminate common mode 
5 signals appearing on the test pair. It is often true that this resistor has a value of 
Z dif /4. The transformer has a limited low frequency response that restricts the use 

of this approach to short loops. 

FIG. 22 illustrates a block diagram of a system which performs broadband 
differential time-domain measurement in accordance with an embodiment of our 

10 invention. FIG. 22 includes waveform generator 2310, sampling head 2510, 
oscilloscope 2320 and loop plant 2100. Sampling head 2510 has input port 2512 
connected to output 2312, loop test points 2514 and 2516 connected to loop access 
terminals 2102 and 2104, respectively, and output port 2518 connected to 
oscilloscope input 2322. Sampling head 2510 includes buffer amplifier 2520, positive 

15 pulse generator 2540, negative pulse generator 2560 and difference amplifier 2580. 
Pulse generators 2540 and 2560 each have an output impedance Z 0 /2 and are 
connected to terminals 2514 and 2516, respectively. The differential impedance 
seen looking into terminals 2514 and 2516 is Z 0 . 

Waveform generator 2310 may, for example be a pulse generator. The 
20 waveform appearing at output 2312 is converted by amplifier 2520 with positive pulse 
generator 2540 and negative pulse generator 2560 into complementary signals which 
are presented to terminals 2514 and 2516. Complementary signals provided at 
terminals 2514 and 2516 preferentially excite a differential mode in loop plant 2100 
and minimize excitation of the several common mode signals. Reflected differential 
25 mode signals appearing at terminals 2514 and 2516 are measured by difference 
amplifer 2580 and absorbed or terminated by the output impedances Z 0 /2. Optimal 
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termination occurs for Z 0 -Z dlff . Difference amplifer 2580 measures the differential 

mode signal appearing across terminals 2514 and 2516 while rejecting common 
mode signals which may be present. Oscilloscope 2320 records the excitation 
waveform and resulting reflected waveforms. Sampling head 2510 accomplishes the 
5 measurement objectives without the inherent low-frequency limitation of the prior art. 
A precisely controlled differential excitation can be applied to the test pairs and the 
resulting echoes can measured with little additional distortion. This greatly simplifies 
the analysis of the measured pair response, as would be done by our method 
described above in processor 2321 . 

10 It should noted that although an oscilloscope 2320 is depicted in FIG. 22 an 

oscilloscope is only illustrative. In accordance with our invention, oscilloscope 2320 
is any device that includes the functions of detecting the signals described 
throughout the present invention, e.g., a receiver, and of displaying such signals. In 
fact, signal display is optional where the operator desires loop make-up identification 

15 or a record. In such an implementation, the display is not necessary and the 
detected signals can be stored in any number of ways, e.g., on a hard drive coupled 
to the receiver, a diskette or CD-rom drive coupled to the receiver, or the capability 
to transmit the data to some other device so that processing can take place. 
Nonetheless, in a TDR type implementation we expect processor means 2321 , the 

20 detection function and the displaying function provided by oscilloscope 2320, and the 
pulse generator 2310 to be included in a single device that may also include our 
sampling head 2510. In addition, it is also quite possible that sampling head 2510 
may be made as a stand alone unit and used as an attachment to present day TDRs. 
FIG. 23 illustrates an implementation of the broadband differential time- 

25 domain sampling head 2510 in accordance with an embodiment of the invention. As 
shown, sampling head 2510 comprises input buffer amplifier 2520, positive pulse 
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generator 2540, negative pulse generator 2560, difference amplifier 2580 and has 
input 2512, positive loop test terminal 2514, negative loop test terminal 2516 and 
measured signal output 2518. Input buffer amplifier 2520 generates complementary 
signals for driving pulse generators 2540 and 2560 and comprises inverter 2612, 
5 level translation 2614 and inverter 2618. Positive pulse amplifier 2540 comprises 
inverter 2620 and series resistors 2626, while negative pulse generator 2560 
comprises inverter 2630 and series resistors 2636. Inverter 2620 has positive bias 
V n 2622, while inverter 2630 has negative V n 2632. 

Inverter 2612 receives input from output 2312 via terminal 2512. The output 
10 of inverter 2612 directly drives the input of inverter 2620 and indirectly drives the 
input of inverter 2618 through level translation 2614. The output of inverter 2618 
drives the input of inverter 2630. The output of inverter 2620 drives loop test terminal 
2514 through series resistor 2626, while inverter 2630 drives loop test terminal 2516 
through series resistor 2636. 
15 Complementary signals appear at the output of inverters 2620 and 2630. The 

outputs of inverters 2620 and 2630 are nominally at zero potential in the absence of 
an input pulse. In one embodiment, input 2512 is provided a fixed amplitude positive 
pulse of duration r , which produces a positive going excursion of duration r at loop 
test terminal 2514 and a negative going excursion of duration r at loop test terminal 
20 P 2516. The complimentary positive and negative going excursions excite a 
differential mode signal on the loop, while minimizing excitation of common mode 
propagation. In one embodiment, inverters 2620 and 2630 comprise high-speed 
CMOS inverters, in which case the positive and negative output pulse amplitudes are 
controlled by the value of the bias voltages V p and V n applied to 2622 and 2632, 

25 respectively. In another embodiment, inverters 2620 and 2630 comprise a plurality of 
individual inverters whose outputs are connected to the respective loop test terminals 
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through a plurality of independent series resistors so as to effect equal current 
sharing among said individual inverters. In yet another embodiment, inverters 2612, 
2618, 2620 and 2630 comprise linear amplifiers with fixed gain such that the 
excitation produced at terminals 2514 and 2516 follow input 2512 in a proportional 
5 fashion. 

The loop signal appearing across the loop test terminals 2514 and 2516 is 
sensed by difference amplifier 2580 comprised of operational amplifiers 2640, 2650 
and 2670, according to known techniques. Operational amplifier 2670 drives output 
2518 through variable resistor R 2680 . Difference amplifier also includes resistors 

10 R 26 42 > R 2652 and variable resistor R 2645 which determine the difference amplifier's 
gain. According to one embodiment, R 2642 and R 2652 have equal values and R 2680 = 
0, in which case the difference amplifier response is 

2R 

V(5 1 8) = ( V(5 14) - V(5 1 6) ) - (1 + 

^2645 

Resistors R 2645 and R 2680 can have fixed preset values. In another embodiment, 
15 the values of resistors R 2645 and R 2680 may be changed at different intervals 

following application of a pulse at input 2512. This permits attenuation of the initial 
excitation pulse and amplification of small reflections appearing later. 

The improvements of the circuitry of the present invention may be better 
appreciated by way of reference to FIG. 24. FIG. 24A depicts an echo response of a 
20 gauge change (3200 feet of AWG 26 cable followed by 3200 feet of AWG 24 cable) 
in accordance with the embodiment of FIG. 22. FIG. 24B depicts the echo response 
of the same gauge change without our common mode rejection circuitry, i.e., using 
the setup of FIG. 20. In FIG. 24B the region noted as A has many peaks, whereas in 
FIG. 24A the same region B is shown as having one peak. Accordingly, our common 
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rejection circuitry acquires a better signal for processing since the additional peaks 
are suggestive of other discontinuities other than a gauge change. 

Accordingly, for the purpose of identifying the make-up of a loop, better 
results will be achieved if the circuitry is used in accordance with the teachings of 
5 section 1 . 

3. Loop Make-Up Identification on a Real Loop 

In this section, an example of our identification method will be given. In 

particular, starting from data measured on an unknown loop, an example on how to 
perform identification following the steps described in the flow-chart of FIG. 17B will 

10 be given. No a priori information on the loop will be assumed, except for some basic 
constraints on the location of discontinuities. These constraints derive from the fact 
that a probing pulse of 5 Volts (differential) is used, and this limits the range of 
identifiable loops. In particular, a set of preliminary experiments allowed us to state 
that, in principle, any topology that satisfies the following constraints can be 

15 unambiguously identified if a 5 Volts (differential) pulse were used: (1) Maximum loop 
length less than 9 kft; (2) Gauge changes located within 5 kft; (3) Bridged taps 
located within 6 kft; 

Obviously, if a stronger pulse were used or some sort of amplification were 
introduced at the receiver these limits would be higher. In particular, it has been 
20 calculated that a minimum amplification of a factor of 50 would be necessary to 
locate accurately any discontinuity located within 1 8 kft. 

Two snapshots of the unknown loop were taken. The first one with a pulse of 
3 /us, and the second one with a pulse of 5 {is. Following the notation introduced in 
section 1, these two snapshots are labeled di(t) and d 2 (0. respectively, and are 
25 plotted in FIG. 25 and FIG. 26. 

The real topology of the loop under investigation is shown in FIG. 27. 
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Following the flow-chart, d^t) is selected as the first data set, or "slice", to 

use. 

Initial Iteration: i = 0 

At step 1708 or iteration ?=0, the gauge of the first loop section has to be detected. 

5 This is done by hypothesizing, step 1709, the four topologies r, (0) , T 2 {0 \ T 3 (0 \ and r 4 (0) , i.e., 
the topology of a loop that is constituted of a single loop section of gauge X (X=26, 24, 22, 
19) and of infinite length. Given the four topologies 7} (0) , (/=1,...,4), the corresponding four 
computer generated waveform h/ 0 \t) (/=1,...,4) are then computed, step 1711. The plot in 
FIG. 28 shows di(t) and the four waveforms hf°\t) (/=1, - • - ,4). 
10 The metric associated with each waveform h/°\t) (/=1,...,4), i.e., the MSE between 

di(t) and the waveforms hf°\t) (/=1,...,4), is now computed as shown in Table 1. Since the 
probing pulse was 3 //s long, the MSE is calculated starting at an instant t s slightly greater 
than 3 //s, e.g. 3.3 jus. Several ending instants t € are considered. 
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Table 1 

15 The MSE pertaining to the waveform h x {Q \t) is always the lowest, so that the chosen, 

step 1712, topology is 7i (0) , i.e. the gauge of the first loop section is an AWG 26. Therefore, 
we pose 7 <0) = T { (0) and h (0 \t) = h^\t). FIG. 29 shows the identified topology 7< 0> after the 
initial iteration, i.e., z=0. 

Next Iteration i =1 

20 The difference e (1) (t) = d } (t) - h (0) (t) is computed, step 1713, and the result is plotted 

in FIG. 30. An echo exists so its rough time of arrival is estimated, step 1714, with any 
method, for example the derivative. The estimate of the time of arrival of the first echo in 
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e (1) (i) is 10.7 jus. The location of the discontinuity that generated that echo is at 3.6 kft from 
the CO as previously discussed. 

The topologies that can be hypothesized, step 1715, depend on the sign of the echo 
and on the previously identified discontinuities. The possible loop topologies look like the 
5 ones in FIG. 31. Since there are several possible topologies, it would be inconvenient to 
calculate all of the possible waveforms and their metrics. So, as suggested in section 1, for 
each family of discontinuities only the two most distant waveforms of the family will be 
computed. The six topologies that are initially considered are shown in FIG. 32. 

FIG. 33 shows the waveform e°\t) and the computer generated waveforms, step 1717 



10 of all six cases shown in FIG. 32. Since the first echo was at 10.7 jus, the MSE between e x (t) 
and [h/ l \t) - h (0 \t)] (/=1,...,6) is calculated, step 1731, starting at that instant. Several ending 
instants t e are considered as indicated in Table 2. 
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Table 2 

The minimum of the MSE is achieved more often and consecutively in the case Al 



15 than in any other case, therefore it is concluded that the real topology of the loop must belong 
to the family (A) shown in FIG. 31, i.e., it is a gauge change. It is worth pointing out that the 
check on the peaks of (Al) and e (]) (t) shows a difference around 20%, a very high value. 
However, in the case of gauge changes, the check on the peak should be performed after the 
fine-tuning, step 1740, of the estimate of the time of arrival as previously mentioned. 

20 The possible topologies for gauge changes are only three so an additional topology is 

hypothesized, a gauge change with a 22 gauge as shown in FIG. 34. The waveform e 0) (t) and 
the computer generated waveforms of all the possible gauge changes cases are shown in FIG 
35. The MSE of the hypothesized topologies (Al), (A2), (A3) of FIG. 32 and FIG. 34 is 
shown below in Table 3. 
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Table 3 

For several time intervals, the minimum MSE is achieved in all three cases so that the 
choice of the right topology might not be immediate. Although in the last intervals the 
minimum MSE is achieved by (Al), in these situations it is better to double check with the 
5 amplitude of the peaks. The percentage differences of the peaks in cases (Al), (A2), and (A3) 
with respect to the peak of e a) (t), are -21%, +105%, and +42%. Certainly, the peak of (Al) is 
the closest to the one of e (1) (t). 

Therefore, the chosen discontinuity, i.e. 7 <1) , is the one labeled Al in FIG. 35; the 
computer generated waveform that corresponds to is h a \t), 

10 Once the topology has been chosen, a fine tuning, step 1740, of the estimate of the 

time of arrival should be performed. This is accomplished by performing a cross- 
correlation between e (1 \t) and the difference [h (1 \t)-h i0 \t)]. A set of cross-correlation 
functions for different time intervals should be considered. This is shown in FIG. 36, An 
analysis of these functions shows that the maximum of the cross -correlation functions occurs 

15 first at +104, then at +166 and finally at +241 lags further from the central lag. The first two 
cross-correlations are symmetric whereas the skewness of the third function is very high, so 
that the value that is chosen is Ai ag =+166. This means that there is an error of approximately 
550 feet in excess in the first estimation of the location of the discontinuity. 

The updated topology 7°* at the end of step i=l is shown in FIG. 37 and the 

20 corresponding waveform is h°\t). The percentage difference between the peak of the 
computer generated waveform h (l) (t) and the peak of e (1) (t) is now -10%. Although this is an 
improvement over the previous value of -21%, there is still a high difference. This suggests 
that another check with the cross-correlation should be performed. The cross-correlation 
between e {l \t) and the difference [h (l \tyh m (t)] (with h (l \t) updated and referring to FIG. 37) 
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is shown in FIG. 38. The time interval over which it has been computed is the same one that 
yielded the value A ]ag =+166, i.e. [10.7 jus, 31 /is]. The skewness of this function is very low 
and its maximum occurs at A Jag =-33. This means that an additional correction of -110 feet is 
necessary. 

5 The updated topology 7 <1) at the end of step i=l is shown in FIG. 39 and the 

corresponding waveform is h (l \t). The percentage difference between the peak of the 
computer generated waveform h (l \t) and the peak of e (1 \t) is now +5.2%. This value may be 
considered satisfactory, so that no more fine tunings of the time of arrival are performed. 

Therefore, the chosen discontinuity at step i=l, i.e. 7 <1) , is the one shown in FIG. 39 
10 and the corresponding computer generated waveform is h (l \t). The waveform h°\t) and the 
data snapshot d } (t) are shown in FIG. 40 for comparison purposes. 

At this point the process proceeds to the next iteration and in accordance with our 
notation /=/+l=2. 

Next Iteration i =2 

15 The difference e (2) (t) = d } (t) - h a> (t) is computed, step 1750, and the result in plotted 

in FIG. 41. The value of e a \t) around 4 l) is less than 10 mV and this confirms that the echo 
was successfully canceled. 

An echo exists in e a \t) so its rough time of arrival is estimated, step 1714, with any 
method, for example the derivative. The estimate of the time of arrival i 2) of the first echo in 
20 e (2 \t) yields ^ 2) =29.2 /is. The discontinuity that generated this echo is located at 9,730 ft from 
the CO or at 6,570 ft from the last identified discontinuity. 

The possible discontinuities that can generate a positive echo in this case are a gauge 
change or the end of the loop. There are only two possible topologies in this case and they 
are shown in FIG. 42. 
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On the basis of the previous considerations, only two waveforms hj \t) (/=1, 2) are 
generated. The plots of e (2 \t) and [/z/ 2) U) - h (l \t)] (/=1, 2) are shown in FIG. 43. Also the MSE 
between e (2 \t) and [hf 2 \t) - ft (1) (0] 0=! 3 ) is calculated and reported in Table 4. 
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Table 4 

5 The MSE for case (B) is always the lowest, therefore it is concluded that T 2 (2 \ i.e. 

case (B) in FIG. 42, is the most likely topology; the corresponding waveform is now 
h {2 \t)=h 2 {2 \t). The percentage difference of the magnitude of the peaks of [h 2 (2 \t) - h (l \t)] and 
e (2 \t) (see FIG. 43) is -26%. This value is somewhat high and, hopefully, the fine-tuning of 
the time of arrival will reduce that value. 

10 The fine tuning of the estimation of if 2) is accomplished by performing a cross- 

correlation between e (2 \t) and the difference h (2d (t)-h (l \t) over several time intervals having 
z* 2) as starting instant. A plot of the cross-correlation functions is shown in FIG. 44. Since the 
cross-correlation functions remain approximately symmetric extending the interval over 
which they are computed, the fine tuning of the time of arrival i 2) can be performed 

15 exploiting the last cross-correlation, the one computed on the interval ^=49 jus]. An 
analysis of this function shows that the maximum occurs not at the central lag but at 27 lags 
further, i.e. there is an error of approximately 90 feet in excess. However, it has to be also 
noted that this new estimate of i 2) may not be accurate because the time interval [0, 50 jus] for 
d\{t) is too short for a complete analysis of an echo starting around 30 jus. This may be 

20 confirmed by comparing the value of e (2) (t) at r=50 jus with its peak; in fact, the value of e a) (t) 
at the end of the observation time is more than 80% of its peak value, confirming that the 
echo did not die out in the observation window. On the basis of the above considerations, it is 
better start the analysis of the second snapshot that has an observation window of [0, 150 jus]. 



-46- 



WO 01/24482 



PCT/US00/27052 



The cross-correlation between e {2 \t) and the difference /z (2) (r)-/z (1) (0 is now 
recomputed using the second data slice d 2 (t) over several time intervals having ^ 2) as the 
starting instant. A plot of the cross-correlation functions is shown in FIG. 45. Since the cross- 
correlation functions remain approximately symmetric extending the interval over which they 
5 are computed, the fine tuning of the time of arrival ^ 2) can be performed exploiting the last 
cross-correlation, the one computed on the interval [if 2 \ ^=89 fjs\. An analysis of this 
function shows that the maximum occurs not at the central lag but at 172 lags further, i.e. 
there is an error of approximately 570 feet in excess. 

The updated topology 7^ at the end of step i=2 is shown in FIG. 46 and the 
10 corresponding waveform is h°\t). The percentage difference between the peak of the 
computer generated waveform h (l) (t) and the peak of e (]) (t) is now +5.2%. This value may be 
considered satisfactory, so that no more fine tunings of the time of arrival are performed. 

At this point the algorithm proceeds to step i=H-l=3, analyzing the data slice d 2 (t). 

Iteration with i - 3 

15 The difference e°\ty=[d 2 (t) - h (2 \t)] is plotted in FIG. 47. The maximum of the signal 

e°\t) in the interval 150 jus] is 10.8 mV, while its energy is 2.9- 10" 6 . These low values 
suggest that no echo is present and that it is necessary to analyze the next data slice. However, 
since the longest data slices was already analyzed it is concluded that the end of the loop was 
reached. 

20 As the above makes clear our method to identify the make-up of loop by processing 

echoes that result from probing the medium with pulses. The above description has been 
presented only to illustrate and describe the invention. It is not intended to be exhaustive or 
to limit the invention to any precise form disclosed. Many modifications and variations are 
possible in light of the above teaching. The applications described were chosen and described 

25 in order to best explain the principles of the invention and its practical application to enable 
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modifications as are suited to the particular use contemplated. 
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CLAIMS 

1 . A system for performing broadband differential time domain measurement 
on a subscriber loop having access terminals comprising: 
5 a waveform generator having an output; 

a detector having an input; and 

a DC coupled sampling head having an input port connected to said 
waveform generator input, two loop test points each connected to the loop access 
terminals, and a output port connected to said detector input. 
10 2. The system of claim 1 wherein said sampling head comprises: 

a buffer amplifier having an input and first and second outputs, said buffer 
amplifier input coupled to said waveform generator output; 

a positive pulse generator coupled to said buffer amplifier first output; 

a negative pulse generator coupled to said buffer amplifier second output; 

15 and 

a difference amplifier having an output and first and second inputs, said 
difference amplifier first input coupled to said buffer amplifier first output and to one of 
the loop access terminals, said difference amplifier second input coupled to said 
buffer amplifier second output and to the other of the loop access terminals, and said 

20 difference amplifier output coupled said detector input. 

3. The system of claim 2 wherein said buffer amplifier comprises a first 
inverter coupled to said buffer amplifier input and said buffer amplifier first output, a 
second inverter coupled to said buffer amplifier second output, and means for level 
translation coupled to said second inverter and a resistor, said resistor being coupled 

25 in series to said first inverter and buffer amplifier first output. 
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4. The system of claim 3 wherein said positive pulse generator comprises an 
inverter in series with a resistor, said positive pulse generator inverter being coupled 
to said buffer amplifier first output, and said positive pulse generator resistor coupled 
to one of said access terminal. 
5 5. The system of claim 4 wherein said negative pulse generator comprises 

an inverter in series with a resistor, said negative pulse generator inverter being 
coupled to said buffer amplifier second output, and said positive pulse generator 
resistor coupled to one of said access terminal. 

6. The system of claim 5 wherein said negative pulse generator inverter and 
10 said positive pulse generator inverter comprises a plurality of individual inverters, 

each having an output that is coupled to the appropriate loop access terminals 
through a plurality of independent series resistors so as to effect equal current 
sharing among individual inverters of said plurality. 

7. The system of claim 6 wherein said inverters comprise high-speed CMOS 
15 inverters. 

8. The system of claim 7 wherein said waveform generator is a pulse 
generator. 

9. The system of claim 7 further comprising means for displaying a 
measurement result. 

20 10. The system of claim 8 wherein said display means and said detector and 

coupled together in an oscilloscope. 

1 1 . The system of claim 9 further comprising a processor having a process for 
identifying loop makeup, said processor being coupled to said waveform generator 
and said detector and said processor tangibly embodying a program, said program 

25 comprising the steps of: 
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transmitting at least one probing one pulse onto the loop from the pulse 
generator; 

acquiring data based on said at least one transmitted pulse at the 
oscilloscope, said received echoes being caused by discontinuities on the loop; 
5 identifying the gauge of a first loop section of the loop based on said acquired 

data; 

calculating a slowly decaying signal based on said identified gauge; 
subtracting said slowly decaying signal from said acquired data to produce 
compensated data, said compensated data having echo signals which represent 
10 discontinuities on the loop; and 

determining the loop make-up based on said compensated data. 

12. A sampling head for performing broadband differential time domain 
measurement on a subscriber loop having first and second access terminals 
comprising: 

15 a buffer amplifier having an input and first and second outputs; 

a positive pulse generator coupled to said buffer amplifier first output; 

a negative pulse generator coupled to said buffer amplifier second output; 

and 

a difference amplifier having an output and first and second inputs, said 
20 difference amplifier first input coupled to said buffer amplifier first output, said 
difference amplifier second input coupled to said buffer amplifier second output. 

13. The sampling head of claim 12 wherein said buffer amplifier comprises a 
first inverter coupled to said buffer amplifier input and said buffer amplifier first output, 
a second inverter coupled to said buffer amplifier second output, and means for level 

25 translation coupled to said second inverter and a resistor, said resistor being coupled 
in series to said first inverter and buffer amplifier first output. 
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14. The sampling head of claim 13 wherein said positive pulse generator 
comprises an inverter in series with a resistor, said positive pulse generator inverter 
being coupled to said buffer amplifier first output. 

15. The sampling head of claim 14 wherein said negative pulse generator 
5 comprises an inverter in series with a resistor, said negative pulse generator inverter 

being coupled to said buffer amplifier second output. 

16. The sampling head of claim 15 wherein said negative pulse generator 
inverter and said positive pulse generator inverter comprises a plurality of individual 
inverters, each having an output that is coupled to the appropriate loop access 

10 terminals through a plurality of independent series resistors so as to effect equal 
current sharing among individual inverters of said plurality. 

17. The sampling head of claim 16 wherein said buffer amplifier input is 
coupled to a pulse generator, said difference amplifier first and second outputs are 
connected to the first and second access terminals respectively, and said difference 

15 amplifier output is connected to an oscilloscope. 



-52- 



WO 01/24482 



PCT/US00/27052 



AMENDED CLAIMS 

[received by the International Bureau on 9 February 2001 (09.02.01); 
original claims 1-17 replaced by new claims 1-16 (4 pages)] 

1 . A system for performing broadband differential time domain measurement 
on a subscriber loop having access terminals, said system comprising: 

a waveform generator having an output; 

a detector having an input; and 

a DC coupled sampling head comprising 

a buffer amplifier having an input and first and second outputs, said 
buffer amplifier input coupled to said waveform generator output; 

a positive pulse generator coupled to said buffer amplifier first output 
and to one of said loop access terminals; 

a negative pulse generator coupled to said buffer amplifier second 
output and to the other of said loop access terminals; and 

a difference amplifier having an output and first and second inputs, 
said difference amplifier first input coupled to said one of said loop access 
terminals, said difference amplifier second input coupled to said other of said 
loop access terminals, and said difference amplifier output coupled to said 
detector input. 

2. The system of claim 1 wherein said buffer amplifier comprises a first 
inverter coupled to said buffer amplifier input and said buffer amplifier first output, a 
second inverter coupled to said buffer amplifier second output, and means for level 
translation coupled to said second inverter and a resistor, said resistor being coupled 
in series to said first inverter and buffer amplifier first output. 

3. The system of claim 2 wherein said positive pulse generator comprises an 
inverter in series with a resistor, said positive pulse generator inverter being coupled 
to said buffer amplifier first output, and said positive pulse generator resistor coupled 
to one of said access terminal. 
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4. The system of claim 3 wherein said negative pulse generator comprises 
an inverter in series with a resistor, said negative pulse generator inverter being 
coupled to said buffer amplifier second output, and said positive pulse generator 
resistor coupled to one of said access terminal. 

5. The system of claim 4 wherein said negative pulse generator inverter and 
said positive pulse generator inverter comprise a plurality of individual inverters, each 
having an output that is coupled to the appropriate loop access terminals through a 
plurality of independent series resistors so as to effect equal current sharing among 
individual inverters of said plurality. 

6. The system of claim 5 wherein said inverters comprise high-speed CMOS 
inverters. 

7. The system of claim 6 wherein said waveform generator is a pulse 
generator. 

8. The system of claim 5 further comprising means for displaying a 
measurement result. 

9. The system of claim 8 wherein said means for displaying and said 
detector are coupled together in an oscilloscope. 

10. The system of claim 8 further comprising a processor having a process for 
identifying loop makeup, said processor being coupled to said waveform generator 
and said detector and said processor tangibly embodying a program, said program 
comprising the steps of: 

transmitting at least one probing pulse onto the loop from the pulse generator; 
acquiring data based on said at least one transmitted pulse at the 
oscilloscope, said received echoes being caused by discontinuities on the loop; 

identifying the gauge of a first loop section of the loop based on said acquired 

data; 
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calculating a slowly decaying signal based on said identified gauge; 

subtracting said slowly decaying signal from said acquired data to produce 
compensated data, said compensated data having echo signals which represent 
discontinuities on the loop; and 

determining the loop make-up based on said compensated data. 

11. A sampling head for performing broadband differential time domain 
measurement on a subscriber loop having first and second access terminals 
comprising: 

a buffer amplifier having an input and first and second outputs; 

a positive pulse generator coupled to said buffer amplifier first output; 

a negative pulse generator coupled to sard buffer amplifier second output; 

and 

a difference amplifier having an output and first and second inputs, said 
difference amplifier first input coupled to said buffer amplifier first output, said 
difference amplifier second input coupled to said buffer amplifier second output. 

12. The sampling head of claim 11 wherein said buffer amplifier comprises a 
first inverter coupled to said buffer amplifier input and said buffer amplifier first output, 
a second inverter coupled to said buffer amplifier second output, and means for level 
translation coupled to said second inverter and a resistor, said resistor being coupled 
in series to said first inverter and buffer amplifier first output 

13. The sampling head of claim 12 wherein said positive pulse generator 
comprises an inverter in series with a resistor, said positive pulse generator inverter 
being coupled to said buffer amplifier first output. 

14. The sampling head of claim 13 wherein said negative pulse generator 
comprises an inverter in series with a resistor, said negative pulse generator inverter 
being coupled to said buffer amplifier second output. 
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15. The sampling head of claim 14 wherein said negative pulse generator 
inverter and said positive pulse generator inverter comprise a plurality of individual 
inverters, each having an output that is coupled to the appropriate loop access 
terminals through a plurality of independent series resistors so as to effect equal 
current sharing among individual inverters of said plurality. 

16. The sampling head of claim 15 wherein said buffer amplifier input is 
coupled to a pulse generator, said difference amplifier first and second outputs are 
connected to the first and second access terminals respectively, and said difference 
amplifier output is connected to an oscilloscope. 
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STATEMENT UNDER ARTICLE 19 (1) 

The new set of claims 1-16 better focus the claims on that which applicant regards as the 
invention. The amendments have no impact on the description and drawing of the application. 
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